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Abstract 

In fluids and plasmas with zonal flow reversed shear, a peculiar kind of transport barrier appears 
in the shearless region, one that is associated with a proper route of transition to chaos. Using 
a symplectic nontwist maps, which model such zonal flows, we investigate these barriers. In 
particular, the standard nontwist map, a paradigm of nontwist systems, is used to analyze the 
parameter dependence of the transport through a broken shearless barrier. Varying a proper control 
parameter, we identify the onset of structures with high stickiness that give rise to an effective 
barrier near the broken shearless curve. Several diagnostics are used: stickiness is measured by 
escape time and transmissivity plots, while barriers that separate different regions of stickiness are 
identified by employing two indicators in the phase space, the finite-time rotation number and the 
finite-time Lyapunov exponent. 
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I. INTRODUCTION 



Internal transport barriers that appear in Hamiltonian dynamical systems have been 
proposed as an explanation for the cessation or reduction of transport in physical systems 
that describe fluids (e.g. jl[ l^j]) and plasmas (e.g. {2]). These barriers may have various 
physical or dynamical origins, yet they can and have been used to control experiments and 
sometimes to improve desired confinement of trajectories. Thus, there is justification for 
studying these barriers in the general context of Hamiltonian systems, which we do here. 

A peculiar kind of transport barrier exists in fluids and plasmas with a nonmonotonic 
equilibrium zonal flow, which give rise to orbit topologies that can only exist with reversed 



shear 
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i.e., with a nonmonotonic rotation number profile. The barriers appear in the 



shearless region of nontwist Hamiltonian dynamical systems and present their own typical 
characteristics with a proper route of transition to chaos They possess robustness - 

persisting even for high amplitude perturbations - and have an effective capacity to reduce 
the transport even after invariant tori are broken 13j. Invariant barriers persist until the 



destruction of the shearless invariant curve 
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14j . but the capacity to reduce transport 



remains and is credited to the stickiness around islands that remain in the shearless region 



Such barriers have been numerically and experimentally identified in several nontwist 
dynamical systems such as those that describe magnetic field lines in toroidal plasma devices 

nn 



with reversed magnetic shear [15|, Il6| , the advection of a passive scalar 



shear flow 



171 ]. travelling waves in geophysical zonal flows 



w an incompressible 
J, the E X B-drift 



motion of charged particles in a magnetized plasma under the action of a time-periodic 
electric field from an electrostatic wave 
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191 ] . and laser-plasma coupling [20] 



The mentioned barrier properties have been theoretically derived for the standard non- 
twist map (SNTM), a paradigmatic example of nontwist systems proposed by del Castillo- 
Negrete and Morrison in 1993 [jjj], and interpreted as a consequence of successive bifur- 



cations of the shearless invariant curve 
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211 ] . This scenario shows for the nontwist 



standard map the relevance of the location of the shearless region where the transport reduc- 
tion occurs. Thus, for all nontwist systems, the transport reduction should be observed in 
the shearless region and not necessarily in high shear regions as for other barriers proposed 
to exist in twist systems 22 1. 
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The standard nontwist map can be regarded as arising from intersections of phase space 
trajectories with a given surface of section. They can also appear as stroboscopic samplings 
of trajectories in a time- dependent system at fixed time intervals. This canonical map is 
convenient for investigating transport, inasmuch we can compute a large number of iter- 
ations in a short time with minimal propagation of numerical error, as occurs with most 
schemes for solving differential equations. This property is particularly important in studies 
of Hamiltonian transport, which require computation of phase space trajectories over very 
long time intervals. 

In nontwist systems, after the shearless curve breakdown, chaotic orbit stickiness is high in 
the shearless region and, consequently, the chaotic transport is reduced in this region. For the 
standard nontwist map, this local transport reduction has been associated with an effective 
transport barrier and characterized in terms of the orbit escape time and transmissivity 28 ]. 



Moreover, it was suggested that the sensitive dependence of these quantities has the same 



parameter dependence as the dominant crossings of stable and unstable manifolds 2 



In this paper, for some parameter ranges for which the transport barrier of the standard 
nontwist map 14j is broken, we identify remaining stickiness structures that reduce the 
transport in the phase space region through the broken shearless barrier. We show how 
these stickiness structures, determined by the homoclinic tangles of the dominant remaining 
dimerized islands, change with the control parameters and modify the observed transport 
within the shearless region. Moreover, we also identify these structures, in phase space, by 
mapping out the finite time Lyapunov exponent (FTLE) and the finite time rotation 
number (FTRN) 9| in the shearless region. The FTLE and FTRN allow us to identify 
structures that separate different regions of stickiness in phase space. Other indicators have 
been used in Hamiltonian systems to delineate regions with other dynamical properties, 
besides the stickiness, as Lagrangian coherent structures |4J and resonant zones 8|, but we 
note the convenience and simplicity of the FTLE for identifying the stickiness structures 
and effective barrier onset. 

The rest of the paper is organized as follows: in Section [TT] we describe two qualitatively 
different transport regimes related to the separatrix reconnection and breakup of shearless 
curves in the SNTM. Section II I II is devoted to a characterization of the transport barriers, 
which are stickiness structures of the SNTM. Section HVl uses the infinite-time rotation num- 
ber to evidence period-three satellite islands just after the formation of transport barriers. 
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Section |V] introduces the finite-time rotation number, an indicator for stickiness structures. 
Section ED uses the finite-time rotation number ridges to visualize transport-related escape 
channels. Finally Section IVIII contains our Conclusions. 

II. STANDARD NONTWIST MAP 

We consider two-dimensional area preserving maps of the following general form: 



where (x n ,y n ) are the normalized angle and action variables, respectively, of a phase space 
trajectory at its nth piercing with a Poincare surface of section. We assume that i/„6l and 
x n G (—1/2, 1/2], and also that f(x n ) is a period-1 function of its argument, representing a 
perturbation strength (whose time-dependence is a periodic delta function). 

The function g(y n +i) is the winding number of the unperturbed phase-space trajectories 
lying on nested tori, its derivative being the so-called shear function. If the function g is 
monotonically increasing or decreasing for all values of interest, the corresponding shear 
does not change sign, which amounts to the following twist condition: 



where c6R. The loci where g'{ys) = 0, i.e. where the shear changes sign, define shearless 
curves in phase space. 

We focus on systems for which the winding number g(y) in a representative map of the 
form ([T]) and ([2]) violates ([3]) by having one extremum at y = y$, at which the shear changes 
sign. In the vicinity of this point the lowest-order approximation for the function g(y) is a 
quadratic one. If, in addition, we keep only one sinusoidal mode in the perturbation term 
of ([1]) and ([2]), we obtain the standard nontwist map (SNTM), introduced in Ref. jjjj : 






(1) 
(2) 




(3) 




(4) 
(5) 




where x E [-1/2, +1/2), y e R, a 6 (0, 1), and b > 0. 
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FIG. 1: (color online) Poincare section of the standard nontwist map and © for b = 0.6 and 
(a) a = 0.364; (b) a = 0.8; (c) a = 0.80552 and (d) a = 0.8063. 

Let us first consider the unperturbed case (6 = 0). The twist condition ([3D is violated 
at the point ys = 0, defining a shearless curve {(x,y)\ — 1/2 < x < 1/2, y = ys = 0}. 
The quadratic form of g around ys = leads to two invariant curves, at y = ±?/o with the 
same winding number a(l — y$) at both sides of the shearless curve. As the perturbation is 
switched on (b ^ 0) two periodic island chains appear at the two invariant curve locations, 
and the former shearless curve becomes a shearless invariant tori separating these two island 
chains. There are also chaotic layers attached to the "separatrices" of both island chains, 
as expected from the presence of homoclinic crossings therein. These chaotic layers are not 
connected, though, as far as there are invariant curves near the shearless invariant tori acting 
as dikes, preventing global transport. 

A representative example is depicted in Fig. (Ha), where a Poincare section of the SNTM 
is depicted for b = 0.6 and a = 0.364 (in the following we shall fix this value of b and vary only 
the parameter a). We observe two island chains with three islands each, with winding number 
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0.3. In the unperturbed map the corresponding invariant curves are located at yo = ±0.42. 
The local maxima of the perturbed winding number profile define a shearless invariant curve, 
whose existence can be inferred between the two island chains. The island chains bordering 
the shearless invariant curve are transport barriers, since chaotic trajectories above and 
below do not mix at all (in Fig. QJa) they have been represented in different colors). 

If the parameters are further modified another noteworthy feature of nontwist maps can 
emerge, depending on the parameter space region. In one scenario (generic reconnection) 
the island chains with the same winding number approach each other and their unstable 
and stable invariant manifolds suffer reconnection. In the region between the chains, there 
appear new invariant tori called meandering curves (which are not KAM tori, though, since 
the latter must be a graph over x, while meanders are not). The periodic orbits remaining 
eventually coalesce and disappear, leaving only meanders and the shearless torus. This 
set is a robust transport barrier, as illustrated in Fig. QJb), where two chaotic orbits on 
different sides of the barrier are kept segregated by a shearless curve [12j . The other possible 



reconnection scenario is nongeneric and invo 
possible in nontwist maps with symmetries 



ves the formation of vortex pairs, which is only 



2l|. 



Further growth of the a-parameter causes the breakup of the transport barrier and the 
consequent mixing of the chaotic orbits formerly segregated on both sides of the shearless 
invariant torus [Figs. UX C ) an d (d)]. However, there are subtle differences between Figures 
[T^c) and (d), which are ultimately related to the invariant manifold of unstable periodic 
orbits (saddles) embedded in the chaotic orbit. 

In order to highlight these differences we resort to numerical diagnostics of the transport 
properties along the globally chaotic layer. The first diagnostic, the escape time, is computed 
as follows: for fixed values of the parameters (a, b) we consider a large number of initial 
conditions in a grid of N P = 2 x 10 6 points regularly spaced in the square [—0.5, 0.5] x 
[—0.9,0.9]. Each initial condition was iterated until the ensuing trajectory crosses either 
one of two boundaries placed at (—0.5 < xb < 0.5, ys = ±2.0). The escape time E is the 
time it takes for a given trajectory to reach any of these boundaries and, since it varies 
widely with the initial condition, we work with an average escape time E. 

Another diagnostic of transport is the transmissivity, which is the fraction of orbits cross- 
ing the (formerly existent) transport barrier. It is computed by randomly choosing a large 
number Np = 4.5 x 10 6 of initial condition on the line (—0.5 < xb < 0.5, y = 1.0) and 
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FIG. 2: (color online) Average escape times E (black full line) and transmissivity T (red dashed 
line) in the standard nontwist map of ([!]) and ([5]) for b = 0.6 and a values in the vicinity of the 
shear less curve breakup. 

iterated the resulting orbits for T = 2.0 x 10 6 times. Then we count the number of orbits 
that cross the broken barrier reaching the y = —1 line. The difference between escape time 
and transmissivity is that the former considers orbits escaping from the region near the 
transport barrier through both sides of it, whereas the latter only takes into account those 
orbits that actually cross the transport barrier. 

In Fig. [2] we depict the average escape time E and the transmissivity T as a sweeps 
through an interval of values in the vicinity of the "last" shearless curve breakup. The 
transmissivity (red dashed line) becomes nonzero only after this shearless curve has been 
broken up [such as in Fig. UX C ) an d (d)], whereas the escape time is finite even before 
this breakup [Fig. Hfb)]. For the parameter value used in the Poincare section depicted 
in Fig. [IJc) the transmissivity is quite high (the large red peak in the middle of Fig. |2]), 
with the escape time being relatively low. By way of contrast, in the seemingly identical 
Fig. [T^d) the transmissivity is dramatically reduced (from ~ 6000 to circa 1000), while the 
escape times remain practically the same. This suggests there is a mechanism preventing 
orbit escape through the former barrier, nevertheless allowing diffusion at both sides of the 
barrier. Hence, even though the last shearless curve has been broken an effective barrier can 
remain, which we explore next. 
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III. TRANSPORT BARRIER 



This effective transport barrier is a consequence of a topological reordering of the invariant 
stable and unstable manifolds of periodic orbits embedded in the chaotic region that follows 
the breakup of the last shearless curve j^]. This chaotic region coexists with the remnants 
of "twin" period-11 island chains [28|. The stable and unstable manifolds of the saddle 
points in the chaotic region therein suffer a reconnection at a value of a between those used 
to obtain Figs. []Jc) and (d). Before this reconnection, the manifolds belonging to the upper 
and lower period-11 island chains cross each other many times (intercrossing), forming an 
escape channel responsible for the high transmissivity displayed by Fig. [U^c), as illustrated 
by Fig. Ela), where the stable manifolds of the upper and lower islands are depicted in red 
and blue, respectively. 

On the other hand, in Fig. [3](b) [which corresponds to the Poincare section of Fig. DJd)] 
the manifolds of the upper (lower) island chain have chiefly crossings with manifolds of the 
upper (lower) chain (intracrossings) , hence diminishing transmissivity, while still allowing 
for diffusion. The changing manifold structure responsible for the local decrease of cross- 
barrier diffusion determines the local stickiness, the characterization of which is the purpose 
of the Sections [TV] and IVT1 

Another illustration of the manifold reconnection forming transmissivity channels is pro- 
vided by the numerical experiment depicted in Fig. HI Here we considered a grid of initial 
conditions and computed for each point the average y- value for a given time t esc = 100. If 
< y > was positive (negative) the corresponding initial condition was plotted in red (blue). 
If there were a perfect transport barrier, like a shearless curve between the upper and lower 
twin chains, then there would be a clearcut separation between points evolving to positive 
large y- values (red) and negative large values of y (blue). After the breakup of the last 
shearless curve, the case of high transmissivity [Fig. Hl^a)] clearly shows the existence of 
incursive fingers of the blue region, showing that there are initial conditions above the up- 
per chain going to negative y through the blue channels. The low-transmissivity [Fig. 11(b)] 
suggests that the manifolds after reconnection act as effective transport barriers, with very 
small diffusion between the colors and few identified incursive fingers. In the next Section 
we consider the FTRN as a means for diagnosing this situation. 



S 




FIG. 3: (color online) Invariant stable manifolds of periodic orbits embedded in the chaotic region 
after the breakup of the last shearless curve for the standard nontwist map with b = 0.6 and (a) 
a = 0.80552 and (b) a = 0.8063. The red and blue curves represent manifolds of the upper and 
lower period- 11 island chain remnants. 

IV. ROTATION NUMBER 

Let x i — y M(x) be a map of the circle S 1 onto itself. If the dynamical system is a 
continuous-time flow, then M can be thought of as a Poincare map obtained through suc- 
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FIG. 4: (color online) Fixed-time average y- values for orbits evolving from a grid of initial conditions 
under the standard nontwist map, with b = 0.6 and (a) a = 0.80552 and (b) a = 0.8063. If < y > 
is positive (negative) the corresponding initial condition is plotted in red (blue). 

cessive intersections of the trajectories with a given surface of section in the phase space. 
The rotation number for the trajectory starting at the point Xq is defined as 



which is lifted to M and IT is a suitable angular projection. According to a theorem of 
Poincare, if / is orientation-preserving this limit exists for every initial condition xq G S 1 
and does not depend on x as well [291 ] . 

As a simple example, let us consider a rigid rotation on the circle S 1 given by M(x) = 
x + w. The corresponding rotation number is u = w. If w is a rational number p/q, where 
p and q are coprime integers, the trajectory represents a period-g orbit of the map M, and 
p is the integer number of times the orbit cycles through the x-direction before returning to 
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FIG. 5: Left: one island of the period-11 primary chain in the chaotic region after the breakup 
of the last shearless curve for the standard nontwist map with b = 0.6 and (a) a = 0.80552; (c) 
a = 0.8060; (e) a = 0.8061; (g) a = 0.8063. Right: corresponding local rotation number profiles at 
x = cross-sections. 

its initial position. If w is irrational, then the ensuing (quasiperiodic) orbit densely covers 
the circle S 1 . The rotation number is not defined for chaotic orbits, for the limit in fl6]) does 
not exist in general. 

The rotation number profiles of the period-11 island chain yield information about the 
topological mechanism underlying the manifold reconnection that we described in Section 
and which creates an effective transport barrier. For this sake we have considered a single 
island of this chain [left panels in Fig. [5] and the local rotation number profiles corresponding 
to cross sections of them taken at x = [right panels in Fig. |5]. Each island of the primary 
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FIG. 6: (color online) Rotation number (in colorscale) as a function of y for x = cross-sections 
and different values of a for the standard nontwist map, keeping b = 0.6. White pixels stand for 
values for which the rotation number is not well-defined, since the corresponding orbit is chaotic. 

period-11 chain of the high transmissivity case [Fig. El^a)] is characterized by having an outer 
secondary period-7 chain and an inner secondary period-3 chain. 

As we approach the point of manifold reconnection [Fig. E^c)] the period-7 chain is en- 
gulfed by the surrounding chaotic sea, and an outer period-10 chain emerges. The inner 
period-3 chain, however, seems to disappear. Its reappearance [Fig. [5](e)] occurs slightly 
after the transport barrier is formed, being also present in the low-transmissivity situation 
[Fig. E](g)]. This transition appears in Fig. [6] marked by the two black dots that indicate the 
period-3 chain death and birth, for a critical parameter a, in different ^/-values. 

The topological changes occurring while the manifolds reconnect can be also appreciated 
from the point of view of the rotation number profiles in the right panels of Fig. or, al- 
ternatively, by the diagram depicted in Fig. [BJ where the values of the rotation number are 
shown in colorscale as a function of y for continuous variation of a-values. As a common 
trend, the formation of transport barrier through manifold reconnection is followed by in- 
creasingly high values of the rotation number. This fact suggests that the transport barrier 
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itself may be somewhat connected with comparatively large values of the rotation number, 
and this suggests a diagnostic based on the rotation number, which we consider in Section 
IVl We also emphasize the appearance of a period-3 satellite island ("gumdrops") just after 
the formation of the effective transport barrier. 

V. INDICATORS 

From the example worked out in Sections HTllIVl we learned there is an effective trans- 
port barrier related to manifold reconnection generating intracrossings, and thus the large 
expansion rates are related to the escape routes avoiding the barrier. From the observations 
of Section IIVI the barrier formation is seen to be accompanied by a localized increase in the 
rotation number - this suggests that the barrier can be related to large values of finite-time 
approximations of the local rotation number. In a previous work [9] we presented the same 
idea for identifying coherent structures. The FTRN is a simpler and computationally faster 
method than the FTLE for one and a half degree-of-freedom systems, because the FTRN 
does not require the evaluation of spatial derivatives or additional differential equations 
and this substantially reduces the computer time. We will pursue both indicators here to 
determine barriers. 

The time-iV finite-time rotation number (FTRN) is thus the time-iV truncation for the 
corresponding iterations of the map M 



In general, ujn, like any truncation, depends on the initial condition. While the infinite- 
time rotation number is not well-defined for chaotic orbits, its finite-time counterpart exists 
for any orbit, chaotic or not. The FTRN measures the average rotation angle swept by a 
trajectory over a time interval T, and thus conveys information about the local behavior of 
trajectories, in the same way as the finite-time Lyapunov exponents does (since the latter is 
the local rate of contraction or expansion). 

In Fig. [7(a) and (b) we show the FTRN and largest FTLE respectively, of a region near 
the period-11 island chain in the high transmissivity case previously shown in the manifold 
diagram of Fig. [3(a). The transport channel provided by the intercrossing of manifolds of 
the upper and lower chains is illustrated by the striations of constant rotation number or 
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FIG. 7: (color online) (a) Finite-time rotation number and (b) finite-time largest Lyapunov ex- 
ponent (in colorscale) as a function of the initial condition for the standard nontwist map, with 
a = 0.80552 and b = 0.6. 

Lyapunov exponent crossing the phase space channels between the islands. 

The FTRN and FTLE plots corresponding to the low-transmissivity case are depicted in 
Fig. |8fa) and (b), respectively. Thus we see that the FTRN is an indicator of effective trans- 
port barriers and stickiness structures that survive in the chaotic region after the invariant 
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FIG. 8: (color online) (a) Finite-time rotation number and (b) finite-time largest Lyapunov ex- 
ponent (in colorscale) as a function of the initial condition for the standard nontwist map, with 
a = 0.8063 and b = 0.6. 

curves are broken. These transport characteristics are essentially related to the stickiness 
and recurrence and have been also described with the FTLE [6J for magnetic reconnection. 
Here we see that the same characteristics can be more easily obtained by calculating the 
FTRN, since comparable results with less computer time give essentially the same picture. 
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Therefore, the FTRN is a "fast indicator" of barriers. 

On the other hand, the FTRN (as well as the FTLE) are not good chaos indicators. In 
the infinite time limit the Lyapunov exponent can be used to distinguish chaotic and regular 
trajectories, but it is not efficient, and the rotation number is defined only for the regular 
(periodic or quasi-periodic) trajectories. However, as seen here a grid of both FTRNs and 
FTLEs in the two-dimensional phase space do reveal structures coincident with regions of 
observed stickiness (around resonances) and effective barriers that are not detected for long 
time observations. We argue in favor of the FTRN, since it is faster with the FTLE, and 
proceed with its further use in the next section. 



VI. RIDGES OF FINITE-TIME ROTATION NUMBER 



FTRNs also provide a convenient way to visualize the escape channels related to inter 
and intracrossing transport. In order to do that, we plot the ridges associated with a field 
of FTRNs, representing crests of local maxima. This is defined as follows: suppose one 
has computed the FTRN field in a two-dimensional region u^{x,y). Then, according to a 
procedure developed by Marsden et al. |30| we superimpose an N x N mesh of equally spaced 
initial conditions (xi, yj), with i,j — 1, 2, . . . N and compute the related Hessian matrix 



H(x,y) 



d 2 u N (x,y) d 2 uj N (x,y) 

dx 1 dydx 

d 2 uj N (x,y) d 2 u) N (x,y) 
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where the corresponding derivatives must be computed for all mesh points, e.g. 
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and so on. The smallest eigenvalue of the Hessian matrix ([8]) is given by 



(8) 



(9) 
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The ridges of the FTRN field are defined as the loci where the following conditions are 
fulfilled: 

Vu N (x,y)-n = 0, A n < 0, (12) 

where 

( dui N (x,y) \ 
aJLl J ' < 13 > 

In Figs. E^a) and (b) we plot the ridges of the FTRN field corresponding to the cases 
of inter and intracrossings, respectively, so as to illustrate the usefulness of plotting the 
ridges for delineating escape channels. The ridges (in blue in Fig. [9]) are plotted with the 
boundaries between positive and negative transport that were previously shown in Fig. H] (in 
red), the latter indicating the escape channels for fast transport. In both cases the ridges 
act as the walls for the escape channels, in such a way that in the intracrossings [Fig. [9(b)] 
the ridges form effective transport barriers, whereas in the intercrossing situation depicted 
in Fig. E(a) the barrier opens into a gateway for transport. Such features are responsible for 
the differences in transmissivity shown in Fig. [2j 



VII. CONCLUSIONS 



One of the distinctive features of nontwist maps, and, in particular, of the standard 
nontwist map, is the capability of developing effective transport barriers, which hamper 
diffusion by means of a trapping mechanism similar to that responsible for stickiness in 
Hamiltonian dynamical systems. These broken barriers are only effective on a timescale of 
the order of the experiment or the observation being conducted, for they are stickiness layers 
of chaotic behavior rather than true invariant tori. That is, one expects some transport to 
occur through these barriers, although on a timescale substantially larger than the typical 
duration of the experiment or numerical simulation. 

Such effective transport barriers can be considered as a type of dynamical structure, 
which are interesting objects to focus attention on for analyzing transport in fluids and 
plasmas, inasmuch they separate regions of qualitatively dynamical behavior. For example, 
when observing the dispersion of pollutants in a marine environment, effective transport 
barriers may separate polluted regions from the surrounding waters on a time scale long 
enough to impede large scale spread of the pollutant. Since such structures change with 
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FIG. 9: (color online) Ridges of the finite-time rotation number field (blue) and invariant manifolds 
(red) for the standard nontwist map, with a = 0.8063 and (a) b = 0.80552; (b) b = 0.80630. 

time, one desires fast indicators that reveal the occurrence, if any, of any internal transport 
barriers. 

For the standard nontwist map we investigated the parameter dependence of the transport 
through the broken shearless barrier. Upon varying a proper control parameter we identified 
the onset of high stickiness structures that give rise to the effective barrier near the broken 
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shearless curve. Once this barrier had been formed, there are two qualitatively different 
phase space regions inside the barrier, which were revealed by by the FTLE and FTRN 
indicators. We also observed the appearance of a period-three satellite island just after the 
formation of the effective transport barrier. 

In the present analysis the computer time for the FTRN is almost an order of magnitude 



less than that for the corresponding FTLE, where we used the well-known method of [31 ]. 
We have shown in this paper the usefulness of both diagnostics on the simplest nontwist 
map presenting such internal transport barriers, but we claim that our results are directly 
applicable to problems of practical interest in geostrophic, atmospheric, and marine flows, 
as well as in the plasma confinement in tokamaks where the plasma has a nonmonotonic 
equilibrium profile. 
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